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ABSTRACT 


Classical plasticity theory for macroscopically isotropic 
materials is briefly outlined. 

An extension of this theory to anisotropic materials, as 
developed by Hill, is discussed. 

The elastic response of hexagonal single crystals is invest- 
igated and a speculative analogy between the yield function and the 
incompressible crystal strain energy function is proposed. Yield and 
plastic deformation under hydrostatic pressure are considered qualita- 
tively for single crystals and polycrystalline aggregates. 

The problem of plastic torsion of a thin walled tube is 
studied and relations are developed for the orientation of the principal 
axes of the deformation and the variation of the anisotropic parameters. 
Results of the analysis are compared with an analysis of the same problem 


by Bailey, et al.* 


*Bailey, J. A.3; Haas, S. L. and K. L. Nawab: Anisotropy in 
Plastic Torsion. Journal of Basic Engineering, March 1972, p. 231. 
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CHAPTER 1 


Classical Plasticity 

Leieeincroduction 

The theory of plasticity, at least for macroscopically iso- 
tropic materials, is a "classical" theory in the same sense as the 
classical theory of elasticity. It has been applied to many problems 
involving metal working processes, limit design, and analysis of stand- 
ard mechanical tests and it has a considerable predictive value. For 
materials which are macroscopically anisotropic or in which anisotropy 
develops with deformation, the theory is less well developed and there 
are relatively few complete solutions of technologically important pro- 
blems. Most solutions which do exist, however, require that the material 
be initially isotropic. Many similarities exist between the theories for 
isotropic and anisotropic materials. Therefore, before beginning a dis- 
cussion of plastic anisotropy it is desirable to discuss the basic con- 


cepts of the classical theory of plasticity for isotropic materials. 


1.2 The Classical Theory of Plasticity for Isotropic Materials 


The classical theory of plasticity was developed from ob- 
servations of the macroscopic behaviour of plastic solids under uniform 
states of combined stress. Since data for metals (and particularly 
cubic metals) have predominated, the general theory is related specifi- 
cally to their properties; although the plastic deformation of other 
materials may be described adequately by the theory. The theory is 


based on the assumptions that: 


i AUTIAND Tie’ 


yvatotzesl4 Isoteesld 

-oet vilsotqoseorsem xot tesol ts ,yitotseslq to qrosds ofT 
eit es 9ensa omsa sit at yrood? "rentoesio" sp et ,eistretam ofqosd 
emsidorq ynam ot betiqqe nsad ead 3] vitottesis to yrosda Isoteasls 
«basta to eteviess bre ,mgtesb jimit ,ssressorg agnbitow [stom gntviovat 
40% .sutev evisotbe1q sldesabtenoo 6 asd 3k brs etast [sokmadoam brs 
yqortosins dotdw nt to oigotjoaias yilssiqoosotosm 978 doldw alstissam 
stod2 bie beqolaveb [low esol at yroors oft ,sokstscroish ditw eqolovab 
~orq tmazroqmt yilsotgolondos) to esotjulos sjalqmoo wet yleviseles 916 
interne edd tad situpess ,zavewod ,Jetxe ob dotdw anoktulos 380M .emeld 
to3 astrosd? siz neswied jatxs astiizalimte yas ,stqortoat yilattiat ed 
-alb & gataniged aroted ,stotatedT .alatisian otqottoabas bas oLqotsjoat 
«foo akaad of? eeunatbh ot oldsiteob at 34 yqotjoains otsesiq to aolesuo 


ssisiretam oiqotzoel sot yttotsesiq to yroeds Isotaesilo afas to atqes. 


elgixejaM ostgotjoal sot yirottaelt to yrosiy Lastaee lo eit Sof 
-do mott begoleveb esw yitottaslq to yxoods Isoteeals oft 
motion rebmy abtioa slzasiq to qwotvaetled atqooeoxcam sf to uantiavies 
yisaluotttag bas) Sitewiict Biab sonké «.esatse beatdmoo to a93s3a 
-ttosqe besefe1 et yrosds Isxsnag od3 , bessakmobertg oved (elstem otdus 
redo 20 wokiamroisb sttaalq od3 dguodtis ;aattseqo7q sheds 03 yliso 
ei yroors odT .yxosdt ods yd yLorsupsbs bodtxoaeb ed yam elebreasm 


:3eH3 anotiqmuses ads oo beaad 


a) the material is a continuum. 

b) creep and thermal phenomena may be neglected (i.e., 
temperatures are assumed to be much lower than the 
melting point of the material). 

c) there is no hysteresis loop in elastic loading and un- 
loading. 

d) there is no size effect. 

One further restriction on the classical theory is that the elastic 
strains are small (of the order of Y/E). This limits the theory to 
moderate hydrostatic stresses. 

The basic observation of an ideal metal undergoing a cycle 

of plastic deformation and loading is as follows. Initially, as the 
body is loaded, it responds elastically with the stress and strain lin- 
early related. Under some critical loading condition the material yields 
or begins to deform plastically, and departs from the original linear 
stress-strain relationship. As the strain increases so too does the 
stress (the material is said to work harden in this case). On unloading 
the material responds elastically and linearly so that the slope of the 
stress-strain curve is the same as in the initial stage of loading. The 
elastic constants are assumed to be unchanged by the deformation as long 
as they are defined with respect to the current plastically deformed state 
of the body. During the deformation, the body undergoes no permanent 
plastic volume change. A further observation is that the yielding is not 
influenced by a moderate hydrostatic stress (i.e., small compared to the 
bulk modulus). 

The mathematical relation which defines the limit of elasticity 


under the system of combined stresses is called the yield criterion, and 
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for the isotropic case it is a function. 


SCS ery = 0 Ci al) 


where Jy2Jo2J3 are the three basic invariants of the stress tensor Oo. 
The requirement that moderate hydrostatic stresses not influence yield- 
ing and the assumption that there is no Bauschinger effect further imply 


that the yield function is of the form 
' ' 


EI, i, 500 | (12) 


i] 


v 
where f is an even function of J, and Jy »Jo 


' 
invariants of the stress tensor Oo: =0 -Itr(o) ibe 


q q 
oJy are the three basic 


There have been many proposals for yield criteria, but the forms which 
are most useful are those due to Tresca and von Mises. 

The Tresca yield criterion states that the material yields 
when the maximum shear stress reaches a critical value. The yield func- 
tion is 

O max-O min = 2k where k = Y/2 
and does not involve the intermediate normal and is unaffected by the 
pressure. The von Mises yield criterion states that yielding occurs 

' ' 


when J, reaches a critical value and is not dependent on Js atak. 


The yield function is 
2 


2 2 2 
£ (614) (01-02) + (02-03) + (03-01) = 6k 
2 


2 2 
+Tyz +Tzx ) 


2 2 2 
(ox-oy) + (oy-0z) + (oz-ox) + 6(Txy 


= 6k2 *(1.2.3) 


where k = Y//Y3 is dependent on the strain history. Physically, the von 
Mises yield criterion can be interpretted to imply that yielding occurs 


when: 


*The subscripts 1, 2, 3, refer to the principal stress system. 
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a) the elastic distortional strain energy per unit volume 
reaches a critical value. 
b) the octahedral shear stress reaches a critical value 


(toct =/2 k). 
3 


The von Mises criterion fits most data better than does the Tresca 
criterion, although, it is more difficult to handle mathematically in 
some problems. The Tresca criterion is often used as a piecewise lin- 
ear approximation to the von Mises yield criterion. 
There is no one-to-one stress-plastic strain relationship 
in plastic deformation; the current stress and plastic strain increment 
are related. It has been deduced {1}* that this relation is the form 
deij? = hog dé (102-4) 


doi; 


i] 
where g and h are scalar functions of the stress invariants J, and J, 


and the strain history. This form of the flow rule predicts that there 
is no plastic volume change and the principal axes of the plastic strain- 
increment tensor and the stress tensor are coaxial. When g = f (or some 
multiple or power of f) certain variational principles and uniqueness 
theorems are valid. These are necessary for the development of the 
mathematical thecry of plasticity. Bishop and Hill {2} have shown that 
this relationship is necessarily true for isotropic aggregates of cubic 


crystals. 


*The numbers in parenthesis refer to references listed in the 
Bibliography at the end of this thesis. 
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1.3 Finite Elastic-Plastic Theory 

Classical plasticity theory as just outlined is restricted 
to moderate hydrostatic stresses. When pressures are large (approach- 
ing several percent of the bulk modulus) there is a finite elastic volume 
change which influences the yielding behaviour of the material. Haddow 
and Hrudey {3} have shown that, if Drucker's hypothesis {4} is valid the 
yield function is of the form 

Gace c)e=.0 Cea al 
where Ts Ge 2y. and 8 = o'/p is the ratio of densities in unstressed and 
stressed states;9 = temperature. Further, if there is no plastic volume 
change, the yield function for isotropic materials is of the form 

5 (Ly 4T, 28) =10 (ies. 2) 


LJ 1 


' ' 
wheret =T -— 1/3tr (T) I and I rt, rl, are the three basic invariants 
! 


of tT. The explicit form of the yield function has not been investigated, 


~ 


but if it maintains the same form as that due to von Mises 


ne ney = 2k* (where the usual summation convention for 
tensors is applied) 
then the von Mises yield function is a limiting case for zero pressure. 
The flow rule is of the form 
DT = ea Cle3y3) 


where Dt is the properly invariant form of the plastic rate of deform- 
- ! 


ation tensor that corresponds to the yield state of stress (g(t ,0) = 0), 


1.4 Conclusion 
This has, of necessity, been a very brief outline of the 
classical theory of plasticity of isotropic materials and is meant 


to deal only with some of the major points of the theory. It has left 
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out development of the variational principles and proof of the unique- 
ness theorem as well as the development of methods of solution of 
special problems in plasticity. Further information can be obtained 


from Hill's book {5} especially the first three chapters. 
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CHAPTER 2 


Plastic Anisotropy 

2.1 Introduction 

Whether as a result of microscopic deformation processes or 
the original structure of the poly crystalline aggregate, anisotropy 
exists or develops in most processes involving plastic deformation 
of metals. The classical theory previously outlined is incapable of 
describing accurately the behaviour in such cases. To construct a 
theory which will deal with plastic anisotropy it is necessary to mod- 
ify the isotropic theory. The anisotropic theory requires a yield 
criterion and flow rule which are directionally dependent and which are 
capable of changing as the state of anisotropy of the element changes. 
One further condition placed on the theory is that for vanishingly small 
anisotropy, the yield criterion and stress-strain increment relation 


should approach those of the isotropic case. 


2.2 Yield Criterion 

The theory as developed here is set up to consider states 
of anisotropy with three mutually orthogonal planes of symmetry at every 
point. The principal axes of anisotropy, defined by the intersections 
of the planes of symmetry, may vary in direction throughout the body 
and are used as local Cartesian axes of reference for elements of the 
body. Following the form given by Green and Naghdi {6} the yield criter- 
ion is assumed to be of the form 

f (1, ekl’) = 0 (27251) 
where° kl and ek1? are the stresses and plastic 


strains respectively referred to some set of axes in the 


body. 
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Although not specifically stated in Green and Naghdi it is assumed here 


that the stresses and plastic strains are referred to the principal 
axes of anisotropy of the element. 

Following the development due to Hill, {7}, the explicit 
form of the yield criterion is assumed to be a quadratic function in 
the stress components which reduces to von Mises yield criterion for 
vanishingly small anisotropy. The yield criterion is of the form 

2£ (oij) = F(oy-oz)* + G(oz-ox)* + H(ox-oy)? 


2Ltyz 2+ 2Mtzx* 


+ 2Ntxy* = 1 G2% 2512) 

where F, G, H, L, M, N are parameters characteristic of the current 
state of anisotropy. 

The yield criterion is this form implies 
a) there is no Bauschinger effect. 
b) moderate superimposed hydrostatic loads do not influence 

yielding. 
If X, Y, Z and R, S, T are the tensile yield stresses and the shear 


yield stresses respectively referred to the principal axes of aniso- 


tropy, then 


1 = Cytol OOF eldest Whang, Zin med 
x? y2 72 x R? 
le H+F, 2462#1+1-1 2M = 1 (Ziaso) 
Y? x2 72 y- 32 
l=F+G, H=1l+1-1, Wel 
ae x ¥2 2 T? 


Clearly only one of F, G, H can be negative and L, M, N are essentially 


positive. For isotropy 


L=M=N = 3F = 3G = 3H (2.2.4) 
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And the yield criterion (2.2.2) reduces to von Mises' criterion with 
2F = 1/y* 

A full description of the state of anisotropy requires a knowledge of 

the orientation of the principal axes of anisotropy as well as the 

six parameters F, G, H, L, M, N (or the yield stresses X, Y, Z, R, S, 


Tr. 


2.3  Stress-Plastic Strain Increment Relation 

More by analogy with the isotropic case than by any physical 
or theoretical justification, a plastic potential flow rule is assumed 
with the potential function and yield function assumed to be equal. 
This leads directly to the stress-strain increment relations, referred 


to the principal axes of anisotropy. 


dex = di. {H(ox-oy) + G(ox-oz)}, deyz = dA Ltyz, 
dey = dd (F(oy -oz) + H(oy-ox)}, dezx = dAMTzx, (2a. 1) 
dez = dA{G(oz-ox) + F(ox-oz)}, dexy = d\Ntxy, 


Since the yield criterion is not pressure dependent, the plastic poten- 
tial flow rule with the plastic potential indentical to the yield func- 
tion predicts no plastic volume change. Note that reversing the signs 
of all of the stresses in a given stress system reverses the direction 
of the plastic strain increments. The principal axes of stress and 
strain-increment are not coaxial unless the principal axes of anisotropy 


and stress coincide. 
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CHAPTER 3 


Some Microscopic Aspects of the Theory 
3.1 Introduction 

So far the theory has been presented from a macroscopic view- 
point with no attempt to correlate macroscopic observations with the 
properties of single crystals. Bishop and Hill {8}, {9} and Hill {10} 
have investigated the plastic deformation of a polycrystalline aggregate 
under the assumption that the only active means of deformation is shear 
along preferred orientations in preferred planes. They showed that for 
cubic crystals (and particularly face-centred cubic crystals) the yield 
function and plastic potential were identical. The same analysis is 
not applicable to non-cubic crystals where additional deformation modes 
are required to produce a general deformation. Several important aspects 
of the theory, however, can be investigated from the basis of the 
elastic response of single non-cubic crystals and their polycrystalline 
aggregates. Hexagonal crystals and the more general trigonal group, 
are the non-cubic crystals which are studied in the following work. 

Under pressure, hexagonal crystals exhibit both volume change 
and deviatoric strain; this behaviour poses the question of whether yield- 
ing and/or plastic deformation are possible under the action of pressure 
alone. There is also the question of yielding of a polycrystalline hexa- 
gonal aggregate under pressure due to the interaction between the in- 
dividual crystals. 

The generalized form of Hooke's Law is 


eij = Sijkl okl (ore, = 35 0) (3.1.1) 
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with the inverse relation 

Oij = Cijkl ekl (or O= Ce) $331. 22) 
where S and C are the fourth order compliance and stiffness tensors. 
The symmetry of the Sijkl and Cijkl in the first two and last two 
pairs of suffixes, as implied by the symmetry of the stress and strain 
tensors, makes it possible to use the matrix notation proposed by Nye 
{11}, summarized by 

tensor notatim ll 22 33 12,21 413,31 23,32 

matrix notation 1 2 3 6 5 4 
This makes it possible to represent the stress and strain systems by 
six-component vectors and the stiffness and compliance tensors by 6 x 6 
matrices. The vectors and matrices, however, are no longer tensors and 
do not obey the tensor transformationlaw. They will be represented in 
this paper by the symbols 

{o] for the stress vector. 

Ce] for the strain vector. 

[sJ for the compliance matrix. 
and Cc] for the stiffness matrix. 


The stress and strain systems in the matrix form are 


{oJ = [oo and Ce] =ffe, 
O02 e2 
03 e3 
Oy ey 
O5 E5 
06 e¢ 


where 01=011 ,02=022,03=03 3 ,04=023505=031 »06=O12 
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The matrix form of Hooke's Law is 


Ce] = (sjJ(cJ 
with the inverse relation 


Col = Cc] Ce] 


3 yar Hexagonal Single Crystals 


Calo) 


CSielne a) 
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The compliance and stiffness matrices for hexagonal crystals 


are 


Cs] a) Sa O12 O13 
Si2 $11 813 
S13 813 S33 


S44 
2(S11-Si2), 


Cc] Cia Cra Cis 
C12 Cia Cig 
Cis) C13 35 
Cyy 
Cyy 


1/2. (CirCi2) 


From Hooke's Law the response to a general stress system is 


fe] = Cs) Co] 


which for the hexagonal single crystal becomes 


Ce] =| Sii Si2 S13 
512511 1's 
Si3 Si3 S33 
S4uy 
Suy 


2(Si1-S12) 


$110, + S1202 + 51303 
S120, + 81192 + $1303 
$1301 + S1302 + S3303 
Sy4y 
S445 
2(S11-Si2) 6 
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For a hydrostatic stress the response is 


Ce] =P 1S11 + Si2 + S73 
S12 % Sit ch Sg 
S13 + Sia + sae 

0 


0 
0 
which has the hydrostatic and deviatoric parts 


(e) = Pi 2541) 025 ers Sr e4 S35 
3 |2Si11 + 2Si2 + 4813+ S33 
2811 + 2Si2 + 4813+ S33 

0 

0 

0 


and 


celea EOed tO) 2 —sOg seo 58 
Spee S1> SS ieee Sets 
—2(Si1 + Si2 - Si3 - S33) 
0 
0 
0 


where bel and le], are the matrix notation forms of the tensors 
1 

tr(e) I and e =e - 1/3tr(e) I respectively. 

For a deviatoric stress the response is 


[Te] =|01 (S11-Si3) + 02(S12-Si3) 
O1 (S12-Si3) + 02(S11-S13) 
Oy" (Sig—59 3)s +: 02S) 3-58) 
SyyOy 
S445 
2(S11-S12)06 


which has the hydrostatic and deviatoric parts 
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14 
and 
Cel = 1]/0) (2811-81 2-281 34533) + 02(-S11+28) 2-281 3+S33) 
3.) G7(-8)1+4281.9-28) 34533). +02 (283 1-S1 2-281,34533) 
O1 (-S11-S12+4S 13-2833) + O2(-S11-Si 2+4S81 3-283 3) 
384404 
384405 
6(S11-S12)06 
Note that both hydrostatic and deviatoric forms of loading produce 
volume change as well as distortion. The values of the coupling con- 
stant (S11+S12-S13-S33) and the volume change constant (2811+2S812+4813 
+S33) are tabulated for a variety of hexagonal and trigonal materials 
in Table 3.1. The values are calculated from data of the elastic 
parameters in the book by Simmons and Wang {12}, except for Tellurium 
which uses data from Schmid and Boas {13}. Although not explicitly 


shown here, the response of trigonal crystals to hydrostatic stress is 


the same as for hexagonal crystals. 


S.3. Stored Energy Function for Hexagonal Crystals 


The stored energy per unit volume, W, is given by 


iis 1{6I" Cs] [0] (3.3.1) 
= 1[o) * [fe] 
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= 1[o1 92 53 Oy Os oe] S1101 + Si202 + $1303 
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TABLE 3.1 


Coupling and Volume Change Constants For Selected Hexagonal and Trigonal 


Materials 
Material S11+S812-S13-S33 25114281 24481 34833 
(10 *ém?/dyne) (10 '*cm?/dyne) 

Be (300°K) 0294 0.2754 
Cd (300°K) -1.3747 2.0734 
aCo -0.0081 0.5271 
Hf =. 0112 0.9206 
Mg -0.0525 a eth 
Re -0.0004 0.2728 
Te 3222 5.14 

Ti -0.0129 0.9324 
OT1 (300°K) 0.2292 2.8103 
Y =0.1145 2.4200 
Zn ei, 2208 1.6921 
tCalcite (Okb) -0.5629 1.3575 
+Calcite (1kb) -0.5644 1.3933 
tCalcite (2kb) -0.5666 eaart 
+Calcite (3kb) -0.5657 1.3790 
+Calcite (4kb) -0.5728 1.3660 
tCalcite (5kb) -0.5774 1.3523 
+Calcite (6kb) -0.5837 1.3364 
Ru -0.0102 ee PPD) 
+Sb -1.3487 2.5814 
+Bi (300°K) =i ined 1.0741 


+Se (300°K) 13.8240 12,8409 


+ Materials have trigonal crystal structure. 


anse. 


LLes . 
Bsts. 
al. 


0 
§ 
1) 
€ 


0 


8 


S+.728 
Ite 


asrt,o- 
8088 .i- 
@63e .0- 
bsd2,0- 
a09¢.0- 
{¢d3¢ .0- 
8ST2.0- 
avyc,0- 
vése.0- 
£0.10:0- 
(BAL. I~ 


Taal . I- 


hE @ Ne: 


eT 

a 

(°00e) Ito 

Y 

as 

(440) ettolsdt 
(dat) ostaisot 
(dAS) satoLedt 
(d#£) s2kofsot 
(dia) sttoladt 
(#2) ottolsot 
(d3i8) extoledt 


we 


.oqisut3a faieyx> Lenogtz? eved aletxszeM t 


16 


1 
W= 1 P*(2S114+281 24481 3+S3 }+2P(-03 ) (S11+S1 2-81 3-Sa3) 
2 


1 


0) ' "5 J t "5 2 
+811(01 “+02 £96. ) +2812(01 O2 -Og “) +251303 


Ug) No "5 
483303 “+Sun(Oy4 “t0s ~) (35362) 


The stored energy function cannot be split into deviatoric and hydro- 
static parts alone unless the mixed term in the expression vanishes. 
This occurs in general only when (S11+S12-S13-S33) = 0, which is the 
condition of no deviatoric strain under pressure or no hydrostatic 
deformation under a deviatoric stress. (The stored energy per unit 


volume for cubic crystals can be split into deviatoric and hydrostatic 


parts). 
In terms of strains, the stored energy per unit volume is 
given by 
=1[e]’ [c] Ce] (G33) 
2 . 
= 1 fei, e2, es, en, 5, ec] Cii Ci2 Cis e) 
2 Ci2 Cii Cis e2 
Cis C13 C33 e3 
Cuy ey 
Cyy es 
1 € 
(C11-Ci2 : 
= liei, G2, 3, €4,e€5, e6,] C11e14C1 2e24€) 3e3 
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aja Incompressible Hexagonal Crystal 


The incompressible hexagonal crystal is a model which leads 
to an interesting result. Incompressibility involves no volume change 
under any combination of stresses. This produces two conditions 


2811 +2Si2 + 4813 +833 = 0 (37451) 
S11 + Si2 — Siz = S33 = 0 (3.442) 


Together, these conditions imply no deformation at all under hydrostatic 
stress and only deviatoric strain under deviatoric stress. Using the two 
conditions it is possible to eliminate two of the five stiffness para- 
meters. Putting S,;3and S33in terms of S,,and Si2, the result is 


Si3 = —(Si1 + Sia) SSB 


$33 2(3811 + S12) (36644) 


Substituting these relations into the strain energy function, it becomes 


W = ae §11(0, 1O5 24803 Vex Ge°) + 281 2(0, gare | oe) 
a 


+ Suu (On7+05°) (3.4.5) 
The strain energy function for the incompressible hexagonal crystal is 
not pressure dependent. It is positive semi-definite in terms of the 
total stress and positive definite in terms of the deviatoric stress. 
The reason for setting up the incompressible hexagonal crystal 


strain energy function is that a hypothetical analogy exists between this 
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form and the yield function. The yield function for isotropic materials 
corresponds to the form of the strain energy function for the incom- 
pressible, isotropic, cubic crystal, and the yield function for aniso- 
tropic materials set up by Hill {18} has the same form as the strain 
energy function for incompressible orthotropic crystals. If this 
analogy is extended to hexagonal crystals, then a yield function for 
materials with hexagonal symmetry would be of the form 


' ' t | ' ! ' 
ZE(Gis) = A (0, <tas, te0, “+20—°) 428 (0) Ge +204 <Ge- ) 


#C (0424052) = 1 (3.4.6) 
In terms of the actual normal stresses, this is 
2£(014) = A(O) “402-4100, -+00102-100103—-100703-606~) 
+2B (=303 “+3010 3+30 203-3010 24306 ~) (3.4.7) 
+3C(o,7+057) = 1 
where the stresses are referred to the principal axes of anisotropy 


of the element with the three direction parallel to the hexagonal axis. 
The above yield function can be formed from the yield func- 
tion 
2£ (014) = F(oy-0z) 2+6(o2z-ox) 24H (ox-oy) ?+2LTyz+2MT2x?+2NTxy2=1 


if the three supplementary conditions 


G=H=A - F ® (5A - 3B)/9 (3.4.8) 
M =N = 3C/2 (3.4.9) 
L= 2F +G (3.4.10) 


are satisfied. These conditions would prevail at the instant of yield- 
ing of a material with hexagonal symmetry, with the x-axis parallel to 
the hexagonal axis, and are comparable to the isotropy conditions. If 


the yield function and potential function are assumed to be identical 
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and the plastic potential flow rule is used, then only certain modes 


of deformation are possible. 


335 Implications of a Pressure Dependent Yield Function 


A major assumption which is used repeatedly in the following 
work is that permanent plastic volume change cannot take place. In the 
range of pressure up to 100 kbar, Bridgeman {14}did not report any 
permanent plastic volume change for any single crystals which he studied 
in his work on the compressibility of solids. Further, the work of 
Alkins {15}, and Hanson and Wheeler {16} on deformation of poly crystal- 
line aggregates showed that volume changes measured were small (<.5%) 
and most likely due to material inhomogeneities such as voids. They 
also measured extremely small volume increases after the initial de- 
crease, which they took to be evidence of expansion of the grain bound- 
aries in response to the mutual constraints imposed by the loading. 

Consider a plastic body, assumed to be an aggregate of rand- 
omly oriented crystals, under the action of a hydrostatic pressure. As 
the pressure increases so too do the internal stresses. For an aggre- 
gate of cubic crystals, the internal stress field is a uniform hydro- 
static pressure since no distortion of the individual crystals occurs 
under hydrostatic stress. For a random aggregate of hexagonal crystals 
however, the internal stress system cannot be hydrostatic since a hydro- 
static stress produces distortion as well as volume change. If each 
grain were subjected to the same hydrostatic stress, the resulting dis- 
tortion would be incompatible. In order that the individual grains fit 
together properly there must be internal shearing stresses in addition 
to the hydrostatic pressure. As the external hydrostatic pressure is 


increased, these shearing stresses eventually will cause plastic deform- 
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ation in some of the individual crystals. The plastic deformation is 
constrained to be of an elastic order of magnitude by the adjoining 
material which is not in a plastic state. Since plastic deformation 

is an irreversible process in the thermodynamic sense, there will be 
energy dissipation during the deformation. Because of this energy 
dissipation, there must be a net volume change in the cycle of applica- 
tion and removal of the pressure. Since there is no plastic volume change 
for single crystals, the volume change which occurs must be elastic. It 
is produced by the residual stresses between the crystals in the aggre- 


gate and is very small. 


3.6 Conclusion 

In work by Schmid and Boas {17} it has been shown that for 
most metal single crystals, including hexagonal crystals, yielding occurs 
when a critical shear stress is reached in the operative glide system. 
For cubic crystals the operative glide system is octahedral and yielding 
can be predicted using the von Mises yield criterion. (Recall that the 
von Mises yield criterion can be interpretted as implying yielding when 
the octahedral shear stress reaches a critical value.) It is purely 
speculation, but it may be possible that the incompressible crystal strain 
energy function used as a yield function implies yielding under critical 
shear in the operative glide system. If this does prove to be the case, 
the strain energy function (2.2.2) for anisotropic materials would corres- 
pond to the most general class of symmetry (orthotropic) and could be 


used safely for all materials no matter what the structure. 
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CHAPTER 4 
Torsion of a Thin-Walled Tube 


4.1 Introduction 
Plastic torsion of a thin-walled tube produces axial strains 
if the tube is unconstrained and significant axial stresses if the 
tube is constrained axially. These effects can be due to the develop- 
ment of anisotropy in the material. Finite plastic strain is con- 
sidered in this chapter; consequently the analysis is based on a rigid- 
plastic model. The wall thickness is assumed sufficiently small com- 
pared to the diameter that the deformation can be taken as uniform. 
Hill {18} and Bailey et al {19} have studied the effects 
of anisotropy using a geometrical representation for the kinematics 
proposed by Hill, Figure 4.9, under the assumptions that: 
(1) the principal axes of the deformation were near those of 
the isotropic case, implying that the effects of the axial 
and tangential strains on the kinematics could be ignored. 
(2) the principal axes of anisotropy and deformation were 


coincident. 


In the following work the kinematics of the problem are 
considered independently of the geometrical representation, using 
the polar decomposition of the deformation gradient to derive an ex- 
pression for the orientation of the principal axes of the deformation 


which includes the effects of the axial and tangential strains. This 
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results in a more rigorous and elegant analysis. 

The deformation is studied for the special cases of free and 
constrained ends and a general set of equations is developed to study 
the variation of the anisotropic parameters with deformation. The 
yield criterion and flow rule proposed by Hill {18} are used and the 
results are compared with an analysis of the same problem by Bailey 


et al {19} using the experimental data from that paper. 


4.2 Kinematics of Deformation 

The torsion of a thin-walled tube can be considered kine- 
matically as an isochoric deformation equivalent to a uniform extension 
followed by simple shear. Referring to figures 4.1 (a), and (b), the 


deformation is defined by the mappings 
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FIGURE 4.1 


A Diagram of the Geometry of Uniform Extension and Simple Shear. 
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23 
where At = lyet,iAa = 1 t+ea,Az = 1 + ez 
The incompressibility condition implies that AaAtAz = 1. Note also 


that the torsion strain is defined on the basis of the deformed length 


of the cylinder. 


4.3 Deformation Gradient 
The deformation gradient F is given by the product of the 
deformation gradients of each of the simple cases. The matrix of com- 


ponents of the deformation gradient is given by 


{r} = {£} *{e} ={9x}{9x) 
~ ff a aX ax 
=|] K 0 At 0 0 
0 i 0 Oke Aa 0 
0 0 1 0 0 Az (453.0) 
=/At KAa 0 
0 Aa 0 
0 0 Az 


4.4 Polar Decomposition of the Deformation Gradient 


The deformation gradient can be decomposed into the product 


of an orthogonal tensor and a symmetric positive definite tensor in two 


ways: 
F=RU (4.4.1) 
where U = (rir)? is the pure stretch followed by a pure rotation. 


F=VR (4.4.2) 
where V = (FF) is the pure stretch which follows a pure rotation. 
The form of the rotation matrix for a rotation about the three axis is 

{R} =|cosa® -sind 0 
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where the angle © is the angle of rotation of the axes measured positive 
in the counterclockwise sense. The rotation matrix is taken to have the 
two dimensional form given, since one principal axis remains perpendic- 
ular to the wall of the tube. The first decomposition is 


{F} = {R} {u} cos Q-sina 0 Wir Uie Urs 
sin a cosa 0 U21 U22 Ui2 


0 0 ni U3, U32 U33 


At KAa 0 
0 a 0 
0 0 Az 


This is a matrix representation for the nine equations 


Ui1 cos® - U2 sina = Xt (4.4.3) 
Ui2 cosa - U2z2 sina = Kha (4.4.4) 
Ui3 cosd — U23 sina = 0 CB Aerie | 
Ui1 sind + U21 cosa = 0 (4.4.6) 
Ui2 sind + U22 cosa = ha (Q04.7) 
Ui3 sina + U23 cosa = 0 (4.4.8) 
U3i = 0 (4.4.9) 
U32 = 0 (4.4.10) 
U33 = Az (4.4.11) 


Since U is symmetric 


Ui3 = U31 = O and Uo23 = U32 = 0 


From (4.4.6) U2i = Uii tan QO. 

Substituting in (4.4.3) gives U,;; = At cos a 

which implies U2, = -At sin a. 

From (4.4.7) Ui2 = (Aa-U22cosa)/sina 

Substituting in (4.4.4) gives U22 = Aa (cosa - Ksina) 
which implies Ui2 = Aa(sina + KeosqQ). 


Since U12 = U21, Aa(sina + Kcosa) = -At sina 
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which implies that tana= -Kia/ (ja + At). (4.4.12) 
Consequently 
U = 1 _ PrtQa +e) wade 0 
VOa +ht)" + (a) ae jatar he OS RE +(Kia) ? 


Similarly, the second polar decomposition gives 


Eat = {v}{R} = 1Vi1 Vi2 Vi3 cosa -sind 0 
V21 V22 V23 sina cosa 0 
V31 V32 V33 0 0 af 


At KAa 0 
0 a 0 
0 0 Az 


which is a matrix representation of the nine equations 


Vii cos® + Vi2 sin® = At (4.4.13) 

-Vi1 sind + Vr2 cosa = Kia (4.4.14) 
Viz = 0 (4.4.15) 
Voi cosa + Voo sina = 0 (4.4.16) 
-V21 sind + Vo2 cosa = da (44.17) 
Ve3 = 0 (4.4.18) 
V31 cosa + V32 sina = 0 (A256 519)) 
-V3, sind + V32 cosa = 0 (4,4.20) 
V33 =Az 


Since V is symmetric 


Vi3 = V31 = O and Vo3 = V32 = 0 


The remainder of the solution is 
Vii = At cosa - Kia sina 
Vio = At sina + KAa cosa 
V2, = Aa sina 
V22 = Aa cosa 


V33 = Az 
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Setting Viz = Vai gives tana = -Kia/ (ia + At) (4.4.22) 
Consequently 
{v} = 1 At (AatAt) + (KAa)? Kia? 0 
V(AatAt)~ + (KAa)~ | KAa ha (AatAt) 0 
: 0 0 Az (AatAt) *+(KAa) * 
{R} = i (AatAt) Kia 0 
V (AatAt)* + (Ka) Ae cera) : 
V(AatAt)~ + (Kia)? 
When Aa = Xz = 1, which is the simple shear case, tana = -K/2 
{R} = 1 2 K 0 
tT 4+K -K 2 ) 
0 0 v4 + xk? 
{u} = 1 2 K 0 
~ V4 + XK? Ka 22k 0 
0 0 v4 + K? 
{v} = i 2+K* K 0 
” V4 + K? K 2 0 
0 0 v4 + K* 
(These relations for simple shear are easily verified). Henceforth, 


let 


B 


Viva + Xt)? (Kaa) 


4.5 Principal Axes of {U} 


The eigenvalues of U are the roots (nN, N2, Ns) of 


de 


At (AatAt) =n KAadt 0 
KAaAt Aa(AatAt) + (Kia) 2-7 0 
0 0 BAz-1) 


mm er 


The third eigenvalue n3 is determined easily by inspection to be 


13 


BAz 


with the associated eigenvector 


COO. ala 


The other two eigenvalues are the roots of 


Caowe) 


nz - {(aatat)? + (Kra) Fn +{radt(AatAt) + Aadt (Kha) f= 0 


(4.5.2) 
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which are 19N25 = { (AatAt) 7+ (Kia) 2+ (Aa*-At*) “+2 (Aa“+At*) (Kia) “+ (Kia }/2 
Let C = V(Ka7-A7t) +2(Aa #7) (Kia) “+ (Ka) 


(4.5.8) 
Then 


TM = (B*4C)/2 and n2(B7-C)/2 
The principal axes of U defined by the eigenvectors are the principal 


axes of stretch before rotation. 


For n Ni, the associated eigenvector is 


(-{Aa*-\t?)+(KrXa)2-C}, 2KAadt, 0) 
For ) = 2, the associated eigenvector is 


(~{Aa*-At?)+(KAa)74+C}, 2KAadt, 0) 


eigenvector for 
DaR02 


eigenvector 
FOr Te =ais 


Figure 4.2 


Orientation of the Principal Axes of {u} 
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Referring to Figure 4.2, 


6, = tn (2) = tan > -(Aa?-At?+(KAa)?-C)/2KAaAt (4.5.4) 
X2 


=f! 
Oa = tan  (\a7-At?+(KAa)74C)/2KAadt (4.5.5) 


i] 
ct 
ie} 
3 

i 
nas 
ais 
wp |e 
Hp 

ee 

I 


tan (0;+02) = tan 6;+ tan 02= C/KAaAt = © 


l-tan6, tan8>, 0 
Therefore 
6, + 65= 1/2 (4.5.6) 
When Aa = At = Az = 1 
ke = paca (-K+V4+K7) /2 (405.7) 
@2 = tan” (K+VEAR) /2 (4.5.8) 


4.6 Principal Axes of V 


The eigenvalues of V are the roots (N1,;N25N3) of 


det | At (AatAt)+(KAa) 7-n Kia’ 0 
Kia? ha (AatAt)-n 0 J=0 (4.6.1) 
0 0 BAz-n 


The third eigenvalue n3 is determined by inspection to be 
N3 = BAz 
with the associated eigenvector 
(OF 07-1) 
The other two eigenvalues are the roots of 
n?—{ (AatAt) *+(Kra) 2}nt{Aadt (AatAt) 7+ (AaAt) (KAa)*} = 0 (4.6.2) 
which are 


niyo = {(AatAt) 24+ (Ka) 24+VCha2=Xt) 242 (02 atd?t) (KAa) 2+KAa) *}/2 


(4.6.3) 


a ete 20.) Ts, 
ea bie 7 Tone . hie oe 
.20) shaisS\ CO" (ataP# Sieh, =! fi \r, 


PE En ree ay x 


== a sta nai — 


(3.2.8) | S\t #99 + »8 | 
peek « th ok aee 

(v2.8) SC ee a) a = 8 

(8.2.3) °\(Saaen) “aaa = <6 


— 


V to asx [sqtontrd 3.4 


90 (effesi 1) e@toox sda ate V to asulsvasgis SAT 


i) ae i=" (phae(sitnh)3h| dob 
(f.d.a) OO 8f O pm (titeh)ek sha 
(fase - 0 0 
ad of not#Seqent yd bentersetab et «| sulsvasgts brhds sfT 2 
shi = eft . 


tojosvasgis betsinoveas sf3 d3iw 
(£ 0 0) | 
to: etoo1 sft 94s eauseealeeh ows talso oelT 
tity, 0 =e? ast CRAM Gaia snattanl ent omen ty | 
sis doksdw 


PRT aR\4S (whay4* (ahtok)} = efheat 


f 
| om _ 


29 


Then 

Ti = (B74) /2 and N2 (B2-C)/2 
The principal axes of V defined by the directions of the associated 
eigenvectors are coincident with the principal axes of the deforma- 
tion. For n = ni, the eigenvector is 

(-{Aa?-At?—-(Kia)?-C}, = 2Ka?, 0) 
For n=n., the eigenvector is 


(-{\a*-\t?—-(KAa) 74+}, 2KAa’, 0) 


eigenvector for 
if) 


eigenvector 
TOR ans 


Figure 4.3 


Orientation of the Principal Axes of V 
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Referring to Figure 4.3. 
1 1 


ti = tan” (=) = tan -(Aa* -At*? -(KAa)?-C)/2KAa* = (4.6.4) 
X2 
fi 
oe 1 
2 = tan (%) = tan (\a?-At?-(K\a) 74+C)/2KAa? (4.6.5) 
X2 
TE 


tan (Y; +2) = @& 


Therefore, ~, + Wo =1/2 (4.6.6) 
When Aa = At = Az= 1 
y= tan” (K+VK?}/2 (4.6.7) 
p2 = See +V44K?}/2 (4.6.8) 


Henceforth, ¢ = 1 in order to coincide with previous work by Hill {20} 


and Bailey et al {21} on the same problem. 


4.7 Relation Between the Eigenvectors of U and V 
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Figure 4.4 


Relation of the Eigenvectors of U and A 
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Referring to Figure 4.4 and the values of the various eigenvectors 


tan (> - 8; ) tan? - tan, 


1 + tand tan§, 


Kha (Aa-At) (Aa?-At?-C)+(AatAt) (KAa)? 
2 (KAa) “Aadt+(Aa “-At--C) --(KAa) *)/4 


= Kaa (Aa-At) (Aa*-At?-C)+(AatAt) (Kia) ? 
4 (KAa) *(AaAt) “+(Aa--At-) --(Aa--At-)C 


= Kia 
AatAt 


= - tana = tan (-a) CHa) 
Therefore, the angle between the eigenvectors of U and Visa, the angle 
of rotation in the polar decomposition of the deformation gradient. 
When Aa = At = Az = 1. 


Tan (¢ -9,) = K = tan (-a) 
2 


4.8 Tensor Transformation Law for Stresses and Strains 
The tensor transformation law in tensor form is 
€ ‘aR = LaiXBjeij,o aR =LaiLBjoij 

where £00 =cos(xa 6 xi) (i.e., the cosine of the angle between the 
new and old directions measuredeve counterclockwise.) 


Expanding these relations for the rotation of axes in Figure 4.5 gives 


dex = sin*} dea + cos*¢ det - 1 dysin 26 (ieee dp 
2 
dey = cos*o dea + sin*o det + 1 dysin 2¢ (4.652) 
2 
dexy = -1 (dea - det) sin 26 + dy cos 29 (4.8.3) 
2 2 


where dy = 2dyat (dy = engineering strain.) and 


ox = sin*d oa + cos“hot -Tat sin 2¢ (4.8.4) 
Oy = cos" Oa + sin*dot +Tat sin 26 (4.8.5) 
Oxy = -1 (oa - ot) sin 2 +Tat cos 2 (4.8.6) 
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The angle $ is measured positively in the clockwise direction. 


a 


$ 


wee 


Figure 4.5 


Orientation of the (x,y) Anisotropic Axes to the (t,a) Major Axes 


4.9 Variation of the Anisotropic Parameters During Deformation 


As the thin-walled tube is plastically strained, the ini- 
tially isotropic material gradually develops anisotropy. The yield 
criterion given by Hill 

2£(oij) = F(oy-oz) *+ G(oz-ox)*+ Hox-oy) 2+ 2Ltyz + 2Mtzx? 

+ 2Ntxy? = 1 
is assumed to hold. In order to describe fully the state of anisotro- 
py in the wall of the tube it is necessary to know the orientation of 
the principal axes of anisotropy and the six anisotropic parameters 
(or the six yield stresses). It is assumed that the principal axes of 
anisotropy coincide at each instant with the principal axes of the de- 
formation as determined by the angle }. 

SinceTyz andtzx are zero for the torsion of a thin-walled 


tube about the y-axis, the parameters lL and M have no effect on the 
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behaviour of the yield function for this problem. The only way that 
L and M can be evaluated is by direct measurement of the shear yield 
stresses. This is very difficult, if not impossible, for the thin 
wall. It is assumed henceforth that the values of L, M and N are 
equal due to latent hardening of the material leaving four parameters 
to evaluate. 
Since 0z = ot = 0 for the loading in the problem, the stress- 


strain increment relations are 


dex = d\{H(ox-oy) + Gox} = sin*¢dea + cos*ddet - dy sin 2¢ 
; (4.9.1) 
dey = dA{H(oy-ox) + Foy} = cos*$dea + sin*ddet + dysin 26 
; (4.9.2) 
dez = dA{Gox + Foy} = -(dex + dey) = -(dea + det) (4.9.3) 
dexy = d\Ntxy ={-(dea-det) sin? + dycos 26}/2 (4.9.4) 


These relations give three independent equations for the four unknowns. 
Since dd = tdy the stress-strain increment relations become 


H(ox-oy) + Gox = (sin*$deatcos*$det-dysin 2o)/tdy (4.9.5) 
2 

H(oy-ox) + Foy = (cos*$deatsin*$dettdysin 26)/tdy (4.9.6) 
2 


Gox + Foy = (dea + det)/tdy (4.9.7) 

Ntxy = {-(dea - det) sin 26 +dYcos 26}/2tdy (4.9.8) 
One further equation can be obtained from the definitions of the aniso- 
tropic parameters in terms of the yield stresses. 

G+H#= 


, Hepa) hres ft (45.959) 
Z2 


1 Al 

x nes 
In order to obtain explicit relations for the variation of 

the anisotropic parameters using the equations developed, it is necess- 


ary to have measurements of the response to the loading. For the tor- 
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Sional loading 9z = Ot = 0. It is necessary to measure two of €a, €t, 
€z as functions of y. It is also required that one of X, Y, Z must be 
known. Usually €a and €t are the easiest strains to measure and Z is 
the only yield stress which can be measured without destroying the 
specimen. (By measuring the hardness of material perpendicular to the 
wall.) It is also necessary to measure the shearing torque and the 
axial load (if the tube is constrained) as functions of y. The system 


of equations used for the remainder of the analysis is 


Gox + H(Ox-oy) = (t0/tdy\'sin2¢deatcosdet-a¥sin29) (4.9.10) 
2 


Foy + Gox = (t*o/tdy) (dea + det) (4.9.11) 
F+G = (to/z)? (4.9.12) 
Ntxy = (t70/2tdy) (-(dea-det) sin 26 +dYcos 26) (4.9.13) 


where F, G, H, N are in the non-dimensional form with TO being equal to 


the yield shear stress at first yield. Also, 


$ = tan ‘1{\t*-\a® + (KAa)* + C} 
2 Kia- 
where C = V(\a--At~)~ + 2(Aa-tAt~) (KAa)~ + (Kial 


If » is assumed not to vary too much from the case where \a = At = Az = 
then 

> = tan 'I{K +744K?} = 7 +1 tan’ K 

Zz 2 


atl 
G2 


The validity of the assumption may be checked by direct comparison 


with full form for > given above. 
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4.10 Tube Not Constrained Axially 
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a) stress system b) orientation of principal axes 


Figure 4.6 
Stress System and Orientation of the Principal Axes of Anisotropy 
For the Unconstrained Tube 
The only non-zero stress is Tat = Tta = T. The stress system 
referred to the principal axes of anisotropy is 
-Ox = Oy = Tsin 26 


T cos 26 (4.10.1) 


TXY 
The yield function (2.2.2) using these stresses becomes 
2£(oij) = Foy* + Gox* + H(ox -oy)* + 2Ntxy2= 1 


<? (sin?26 (F + G + 4H - 2N))=1 


implying 
t= (N+ (F +6 + 4H - 2N) sin? 26)7? (4.10.2) 


The system of equations for the unconstrained tube is 


G + 2H = - (to/T)*(sin?2p dea + cos* dea - 1 sin2)/sin2¢ 
ss (4.10.3) 
F - G = (to/t)?(dea/dy + det/dy)/sin 2¢ (4.10.4) 
F +G = (to/z)? Ce 10s) 
N = (to/t)” {-(dea - det) tan 2 + 1}/2 (4.10.6) 


dy dy 
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Since €a and €t are functions of Y only, the ratios dea/dy and det/dy 
are equivalent to the derivatives. From the curves published by Bailey 
et al {21} the strains are given by the approximate functions 


€a(y) = (.01) (.25y + .171y? + .025y° + .00104y") (4.10.7) 


tt 


EEC) = -ea(y)i/3 
and the derivatives are given by 


dea(y) = (.01) (.25 + .342) al. O75y> +6200426y°%)) veGe10.8) 
dy 


det(y) = - dea /3 
dy dy 


These functions are fifth order polynomial expansions (note fifth order 
term = 0) and are not valid for y>10 where the functions begin to diverge. 
The expressions for T/ToO and z/Tto proved to be difficult to model and 
the calculation of the variation of the anisotropic parameters was under- 
taken using point by point input of values derived from the paper by 


Bailey et al {21}. The results are presented in Figure 4.7. 


4.11 Tube Constrained Axially 


Though this case was not analyzed or experimented on by 
Bailey et al {21}, it was decided to include it for completeness. The 


stress system and orientation of the principal axes are shown in Figure 


4.8. ne 
v c z 
Oa 
a) stress system b) Orientation of principal axes 


Figure 4.8 


Stress System and Orientation of Principal Axes of Anisotropy for the 
| Constrained Tube 
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Figure 4.7 


Variation of anisotropic 
parameters in torsion of 
a_thin-walled tube, 


—.— Values from paper by Bailey 
et al. 


— — Calculated values from yield 
and shear stresses in paper 
by Bailey. 
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The stress system referred to the principal axes of anisotropy is given 
by 


Ox 


i] 


T(A sin? - sin 2¢) 


T(A cos*} + sin 20) (ata) 


Oy 
Oxy = T(-A sin 26 + cos 26) 
Z 
where A = Oa/T 


Under this stress system the yield criterion (2.2.2) becomes 


.2£(01j) = Foy* + Gox? + H(ox-oy)? + 2Ntx*y = 1 


F(A cos*o+sin2¢) 7+G(Asin*¢-sin2) *+H (-Acos2-2sin2¢) ” 


+ 2N (-A sin 26 + cos 26)* = 1/1? (4.11.2) 
2 


Two special cases were considered 
A) dea = det = dez = 0 


This gives 


= 
ay 
ct 
i) 
=] 
“A 


2 Pe 


and the system of equations becomes 


-G (A sin*$-sin2o) + H(A cos2¢+ 2sin2) (to/t)*sin 26/2 


(ate3) 

F (A cos*$+sin 26) + G (A sin*$- sin2¢) = 0 (Canes) 

FE + G = (to/z)? (V1.5) 

N = (to/t)* cos 26 (Aad U76) 
co 2o+cos 26) 


B) dea = 0, det = -dez 


The stress system remains the same as in the above case. The 


angle » is given by 


@ = tan + 1{(At?-1) +K? + V(d-\t 2) “+2 (14At) 7K74K" } 
2 K et ee. 
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The system of equations is given by 


-G(A sin*-sin 2) + H(A cos 2+ 2 sin 26) 


= uP fede 2 o+ sor gies) gee ay 
¥ 


F(A cos’ + sin 26) + G(A sin*$-sin 2¢) 


= -(t0/T)* dez (4.11.8) 
dy 
F + G = (to)? (4.11.9) 
Zz 
N = (to)* -sin 2odez + cos 2¢ (4.11.10) 
E d 


-Asin 26 + 2cos 26 

There should be no difficulty involved in actually running an experi- 
ment to determine the variation of the anisotropic parameters. The 
only measurement not required before is one of axial thrust produced 
during the torsion. It would also be possible to do a complete analy- 
sis for arbitrary axial stress (as long as it was below the yield point) 
and torsion. (in effect Hohenemsers problem of combined tension and 
torsion could be analyzed for the variation of the anisotropic para- 


meters ) 


4.12 Conclusions 

The paper by Bailey et al {21} which has been cited often 
here had much the right approach to the problem, but the analysis is 
incorrect. The error arises in the transformation equations used: 


dex = 1(dea + det) + 1(dea -det) cos 28 + 1 dysin 20 
v4 2 2 


dea cos?@ + det sin*0 + dy sin 20 (Galedh) 
2 


dey = l(dea + det) - 1(dea-det) cos 20 - 1 dy sin 20 
2 2 2 


2 
dea sin?@ + det cos*6 - dy sin 20 Cato) 
2 
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Hats 


(8.11.4) gab *¢rN6r)- = 
Yo 
(@. 11.4) * (on) =o+1 
s 
(OL. £1.) $$ aon + saboS ma- * (Ox) =v 
b » 


of aoss + OS AleaA- 


etraqxs ne gains yilevtos at beviovnt eaupkvirb on sd bivode ‘sted? 
sfiT .exstematsq otgoxtoetns sda io aolsatiav of? sntims2eb of 329m 
 beowbotq teurls [skxs to sno st aroted bettupss ton temo Tyassm xine 
~yiens sielqmon 5 ob of afdtaeog ad opts bivew 41 .aoterto? oA& gatveb 
(antoq blety sis woled asw 3t ee anol as) sesria isixs ysetttdxs 102 e¢e 
bre nokened bentdmos to meldorq etoamenadoH Joetis aI) .morexro3 bas 
~steq otqotsoaitns sit to motsetiav ond xoi basyfsns ed bluos molete3 


Letetem 


anotauloxo9 SL. 
feito betts mead esd dotdw iS} io 45 yeliad yd 1sqsq oT 
et eteylsne edt tud ,meldorq oft o3 dosomqgs sigis sd3 ihosme bed oxed 
thseu enokdevpo solteorotenss3 sii nt esetis roxxe adT . Soesvtoanl 
OS atayh 2+ 86 dod (a9b- dab). + (396 + s9b)f = xab 
(L.SL.a) OS ake tb + Onnte gab + OS e02 edb = 
0S mite v8 E - 8 woo (a9b-n9b)L ~ (aah OBB) = vb o 
| | s 
(S.SE-8) Saiepee- Tage sabi ate ash = 


@ 


40 


dexy = - 1(d€a - det) sin 28 +:dy cos 28 (4, 12.3) 
2 2 


where "0 is the angie of inclination of the (x, y) anisotropic axes to 

the (a, t) major axes, and is considered positive anticlock- 

wise." (see Figure 4:5). 

Although not explicitly stated, Bailey uses 9 = -¢. The 
correct forms of the transformation equations are (4.8.1), (4.8.2), 
(4.8.3) (they can be verified by using the relations in appendix IV of 
Hill {22}). The error is due to the fact that the (a, t) axes are not 
a right-handed system of co-ordinates, and the correct forms can be 
produced if the a and t are interchanged in the transformation equations 
used by Bailey. 

Also, the equations for ¢ and dd/dy are unnecessarily com- 
plicated; for the isotropic case 


b= sin ? sin p/2 (45, 12.54) 
(1 + (2 cosW siny/2) - 2 cost) 


where tan W = 2/y = 2/K 


dd = (-2  ) E (a-sin'¥/2)™) x (A 12 35) 


dy +4 Z 


x E cosw/2 - A? cinv/a f asta /2(centcost2 - 2siny/2 sec“) 


tan tan“ 
+ 2 siny f ] 


Where A = {1 + (2 siny/2)* - 2 cos} 
tan 


whereas the analysis using the polar decomposition gives (for the 


isotropic case) 


$= tan! 1(y + Waty”) =1 + 1 tan 'y (4.12.6) 
2} hot 2 2 
and 
do= 1 (4.12.7) 
dy Hh 
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(It is possible to get the last relation for » directly from the geo- 


metry of Figure 4.9.) 


| Axis of Cylinder 


Figure 4.9 


Geometry of Deformation for the Plastic Torsion of a Thin Walled Tube 


Another fairly minor point is that it is totally unnecessary, 
and perhaps even harmful in a numerical sense, to put €a and €t as func- 
tions of >. It is more natural to derive the strain functions in terms 
of y and it is imperative that they be in that form if the correct form 


for $ is to be used 
> = tan 11{\t?~-ha? +(yAa) 2+ V(Aa7-)7t) 7+ 2 (ha7+At”) (yAa) *+(yAa) * 
2 via 


For the case studied, the difference in the angles ¢ between the assumed 
fsotropic form and the correct form above is the order of hae radians 
and decreases for very large values of y , showing that the assumption 
is valid. Referring toFig. 4.7, the variation of the anisotropic para- 
meters as calculated by the two analyses is presented. Differences as 


large as 40% in the value of F (aty = -5) and 15% in the value of N 
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(between Y = 1.5 and 2.5) exist between the calculated values, although 


at large strains the differences are less than 5%. The differences are 


most probably due to the use of the incorrect transformation equations 
by Bailey, {19} , although round off and other numerical errors from 


the clumsy forms for > , dd and d€a and det may also have an effect. 
dy dy dy 
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